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detail, in particular how to impose cuts to maximize signal-to-background ratio. 
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I. INTRODUCTION 



Exclusive double diffractive production (EDD) of the Higgs boson has been suggested 
some time ago as alternative to inclusive measurements [j], Hf- Exclusive diffractive dijets 
production attracted recently a lot of attention due to new data from CDF run II Q. The 
standard approach for the calculation of central dijets production in proton- (anti)proton 
collisions is based on the Kaidalov-Khoze-Martin-Ryskin (KKMR) QCD mechanism which 
was initially developed for the central exclusive Higgs production in Refs. Q which is ex- 
pected to provide a really robust signal due to a clean environment and highly suppressed 
backgrounds (see, e.g. Refs. j5|-0|)- For more details on the central exclusive processes and 
related physics, we refer to the most recent reviews in Ref. 0. 

It is known, however, that the process pp — > pp{qq) is dominated by the non-perturbative 
region of gluon transverse momenta, and even perturbative ingredients like the Sudakov 
form factor are not under full theoretical control [||. Uncertainties on exclusive diffractive 
production of Higgs at the LHC were discussed in Ref. [l(| together with uncertainties on 
gluonic jet production which was measured by the CDF collaboration. The problem becomes 
even more pronounced when considering the irreducible backgrounds in central exclusive 
production of Higgs boson originating from the direct exclusive bb pair production in a fusion 
of two off-shell gluons. In particular, in Ref. it was shown that the central exclusive 
production (CEP) of bb jets at LHC (see Fig. [I]), may totally shadow the corresponding signal 
of the Higgs boson in the bb channel (see Fig. [2]), which may lead to significant problems in 
experimental identification. 

Therefore, it becomes very important to investigate the exclusive quark jets production 
in different kinematical domains and quantify the related theoretical uncertainties. On 
the other hand, the analysis of various differential distributions and experimental cuts in 
considered four-body reaction pp — >■ p + "gap" + (qq) + "gap" + p could help in a reduction 
of the corresponding backgrounds. 




FIG. 1: Direct central exclusive bb 
pair production in /cj_-factorization ap- 
proach. It is considered to be the main 
irreducible background for Higgs CEP. 



FIG. 2: Central exclusive production of 
Higgs boson production with its subsequent 
decay into bb pair, which competes with the 
direct bb pair production shown in Fig. [TJ 



Unpolarized exclusive c and c jets production was investigated numerically in Ref. [12]. 
It was found that the whole process is dominated by quark/ ant iquark production with 
low transverse momentum k± in the very forward limit of outgoing protons, whereas it is 
strongly suppressed at high fcj_'s, much stronger than in inclusive case. The same should 
hold for b and b jets CEP, which constitutes the major part of the irreducible background 
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for central exclusive Higgs production 1 . As was demonstrated numerically in Ref. such 
a background turned out to be dominated by small gluon transverse momenta ~ q± (which 
are the same for both active and screening gluons in the forward limit) coming into hard 
subprocess amplitude g*g* — > qq contracted with the gluon transverse polarization vectors 
~ g^/xy'i and integrated over in the diffractive amplitude. The presence of the screening 
gluon in the loop, actually, violates the well known J z = selection rule jl4f . which was 
initially established in one-step mechanisms, like 7*7* and PP fusion processes. Such a 
violation manifests itself in the fact that the hard subprocess amplitude is non-zeroth and 
proportional to gluon transverse momentum, and it is indeed strongly suppressed only in 
the high quark transverse momentum limit, whereas in the low-k± limit it may lead to a 



significant contribution [11, |12 



These observations still suffer from a lack of solid theoretical background. In order to 
predict observable signal from quark dijets at LHC and to make a decisive conclusion about 
irreducible background for Higgs CEP, it is worth to analyze carefully different kinematical 
limits both numerically and analytically. The main goal of this paper is to derive explicitly 
the hard subprocess amplitude g*g* — » qq for any quark helicity states X q and X q in diffractive 
kinematics in two different configurations of large k± 3> m q and small k± <C m q quark 
transverse momenta focusing primarily on large invariant mass of qq dijets. This would give 
us an opportunity to analyze different kinematical asymptotics of the diffractive amplitude. 
Having such amplitudes we can then numerically evaluate differential distributions in quark 
transverse momenta, rapidity, relative angle between the quark jets and invariant mass of 
the qq dijet. These theoretical elements are necessary for upcoming Higgs searches and 
diffractive dijets measurements at LHC. 

In the present paper we extend earlier studies related to Higgs and jet production j3.@.[l~5| 
to the production of quark- ant iquark jets. In our approach we use unintegrated gluon 
distributions as proposed by the Durham group Q. Slightly different gluon distributions, 
fitted to the HERA data, have been used in Ref. [||. The choice of gluon distributions 
brings uncertainties of a factor of about 2. In our case we consistenly use the same gluon 
distributions for the signal and background. In the case of Higgs production we shall use 
off-shell matrix element for g*g* — > H compared to the on-shell matrix element used before. 

The paper is organized as follows. In Section II we consider the general kinematics of the 
central exclusive dijet production. Section III is devoted to a discussion of the diffractive 
amplitude, in particular, its hard and soft constituents. Explicit derivation of the helicity 
amplitudes for the hard subprocess part g*g* — > qq in general kinematics and in some 
important limits is given in Section IV. Electromagnetic 7*7* contribution is discussed in 
Section V. In Section VI we reexamine the central Higgs production taking into account 
gluon virtualities and explore contributions of the exclusive bb pair and Z production as 
backgrounds for Higgs. Section VII contains discussion of numerical results. Finally, some 
concluding remarks and outlook are given in Section VIII. 

II. KINEMATICS OF THE CENTRAL EXCLUSIVE DIJET PRODUCTION 

Inclusive heavy quark/antiquark pair production in the framework of the /c_i_-factorization 
approach [l6[ was considered in detail in Refs. [l7-19|. In particular, it was shown 



1 Other backgrounds although in principle irreducible can be large 13 1 
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that the combination of the /cj_-factorization approach and the next-to-leading-logarithmic- 
approximation (NLLA) BFKL vertex in Quasi-multi-Regge kinematics (QMRK) [20] to- 
gether with the concept of unintegrated gluon distribution functions (UGDFs) gives quite 
good agreement with data on inclusive heavy qq pair production. 

It looks quite natural to apply similar ideas to exclusive diffractive qq production in 
proton- (anti)proton collisions at different energies. In Figs. [3] and H] we show the general 
kinematics for the process pp — > p + "gap" + (qq) + "gap" + p under consideration at the 
parton and the hadron levels, respectively. 



'In o 



hi 




ff(x',x 2 ,qt ± ,qj ± ,i4;t) 



FIG. 3: General kinematics of exclusive diffrac- FIG. 4: Cross section of 2 -)• 4 pro- 
tive qq pair production in pp collisions at the cess of exclusive diffractive qq pair 
parton level. production at the hadron level. 

The decomposition of gluon momenta into longitudinal and transverse parts in the high 
energy limit in the c.m.s. frame is 



Qi = XiPi + Qi±, Q2 = x 2 p 2 + <?2±, < x 1<2 < 1, 

go = x'lPl + X2P2 + Q0±, x'l ~ x' 2 = X < 2i )2 , 9o,l,2 - 90/l/2±- 

Making use of conservation laws 

qi = Pi - p'i - qo, q2=P2-p' 2 + qo, qi + q2 = h + k 2 , 

we write 



s X\X 2 = + |k_L 



Mi-,, Mi- = (k x + h" 2 



M 



qq 



2) , 



Xl,2 



(2.1) 
(2.2) 

(2.3) 



where M q q and y q q is the invariant mass and rapidity of the qq pair, respectively, and 

k± = -(p'i± + p' 2± ) = qi± + q2±. = hx + fa±, 

is its transverse momentum, where q±/2± and ki/ 2 ± are gluon and quark transverse momenta 
with respect to the c.m.s. beam axis. In analogy with Eq. (12. ip . we can write 

P'i = C1P1 + Pi±, P2 = C2P2 + P 2± , Ci,2 = 1 - ^1,2 , (2.4) 

where p'1/ 2 ± = ti,2 in terms of the momentum transfers along the proton lines ti >2 . 

In the c.m.s. frame it is convenient to choose the basis with z-axis collinear to the proton 
beam. Then the proton momenta are 



Pi 



1, 0, 0, 1), 



p 2 =^(l, 0, 0, -1). 



(2.5) 
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Let us choose the y-axis in such a way that qf = —q\ = q y . In these coordinates the gluon 
transverse momenta are 

Qix = (0, gf, q y , 0), g 2± = (0, £, -q\ 0) . (2.6) 

Conservation laws provide us with the following relations between components of gluon 
transverse momenta and covariant scalar products 

, g|i + (gi±g 2 ±) q%± + (gi±g2±) „ tf _ ViL^Il ~ (H±Q2±) 2 

Qi- ^ , 92- jj^j , 9- jj^j sign(g), 

k l = -I k ±| 2 = Qi± + Q2± +Hqi±Q2±), qf/ 2 ± = -|qi/2±| 2 , (2.7) 

where |k_J is the gg-pair transverse momentum with respect to z-axis. On the other hand, 
momentum conservation p' 1± + p' 2± = k± leads to a useful relation 



-ti -t 2 + 2 v / M 2 cos$ = |kj , (2.8) 

where $ is the relative angle between the outgoing protons. 

The appearance of the factor sign(g s/ ) guarantees the applicability of Eq. (12. 7p for both 
positive and negative q y . Note that under permutations gu_ <-> g 2 j_ implied by the Bose 
statistics the components interchange as qf <-> gf and g y <-> — g y . 

In analogy to Eq. (12 .ip one can introduce the Sudakov expansions for quark momenta as 

ki = x\p\ + x q 2 p 2 + h±, k 2 = x\pi + x q 2 p 2 + k 2 ± (2.9) 

leading to 

x i,2 = x\,2 + x \,2 = ^7^e ±yi , 2 = —^e ±y2 , m\ /2L = m\ + |ki/ 2 _L | 2 , (2.10) 

in terms of quark/ ant iquark rapidities j/i, y 2 and transverse masses mu_, m 2 j_. In the 
considered coordinates we write in analogy to Eq. (I2.6P 

ki± = (0, kf, k y , 0), k 2L = (0, kl -k y , 0), (2.11) 

with components satisfying the relation kf + k 2 — qf + gf- By construction, in order to get 
the diffractive amplitude in the covariant form useful in any coordinates, we should relate 
components kf, 2 and k y with the scalar products in the similar way as for gluon momentum 
components (see Eq. ( 12. TP ) : 

fc x = kf ± + (k 1± k 2L ) ^ ^ = k 2 2± + (k 1± k 2± ) ^ ^ = ^/k 2 1± k 2 2± - (fc 1± fc 2± ) 2 g . gn(F) ^ 

(2.12) 

In subsequent calculations we will construct the gg diffractive amplitude in explicitly 
covariant form and analyze its behavior in different regions of the 4-particle phase space. 
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III. DIFFRACTIVE AMPLITUDE 



Generally, in the case of the central exclusive production (CEP) with the leading pro- 
tons, the central system X should necessarily be produced in the color singlet state, such 
that the proton remnants and the X system are disconnected in the color space and their 
hadronisation occurs independently giving rise to rapidity gaps j2lj. So, without the loss of 
generality we are concentrated on the simplest case of qq pair produced in the color singlet 
state. 

According to the KKMR approach 0-0] we write the amplitude of the exclusive diffractive 
qq pair production pp — > p(qq)p as 

where the transverse momenta and the longitudinal fractions of gluons are defined in the pre- 
vious Section, X q , Xq are the helicities of heavy q and q, respectively, f° G is the unintegrated 
gluon density function (UGDF) and V\ q \^ is the hard subprocess g*g* — > bb amplitude. 
Averaging over color indices c±, c 2 of t-channel fusing gluons is made explicitly. The nor- 
malization convention of this amplitude differs from the KKMR one by a factor s. The 
amplitude is averaged over the color indices and over the two transverse polarizations of the 
incoming gluons. The bare amplitude above is subjected to absorption corrections which 
depend on collision energy and typical proton transverse momenta. We shall discuss this 
issue shortly when presenting our results. 

A. Matrix element of the hard subprocess g*g* — > QQ 

Let us consider the subprocess amplitude for the qq pair production via off-shell gluon- 
gluon fusion. The vertex factor V^ 2 = Vj^^ki, fe) in expression (13. L p is the production 
amplitude of a pair of massive quark q and antiquark q with helicities X q , Xq- and momenta 
ki, k 2 , respectively. Within the QMRK approach |20| we have 

TOfli, ft) = nfaKS^' n T = ( 3 - 2 ) 

*- J p,cms 

i,k 

where E vcms = yfs/2 is the c.m.s. proton energy, t c are the color group generators in 
the fundamental representation, u(kx) and v{k 2 ) are on-shell quark and antiquark spinors, 
respectively, b 11 ", b^ u are the effective vertices (13. 3p arising from the Feynman rules in the 
QMRK approach illustrated in Fig. [5j 

y (h , ki) = m l, r - , (3-3) 



(qi - h) 2 - m? (h + k 



2) 



b^(k 2 , kl ) = 7 ^ 1 ~k + "V - ^ 



(qi - h) 2 - ml (ki + k 2 ) 2 

where by the Dirac convention a = 7 • a for any 4-vector a M is adopted, T^ u "{qi,q 2 ) is the 
effective three-gluon vertex. These effective vertices were initially proposed for massless 
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quarks in Refs. [20 



The effective ggg- 



and then extended for massive case in Ref. (13, 
vertices are canceled out when projecting the qq production amplitude Eq. (13. 2 j) onto the 
color singlet state, so only the first two diagrams in Fig. [5] contribute to the final result for 
the production amplitude. Since we will adopt the definition of gluon polarization vectors 
proportional to transverse momenta qy 2 ±, i-e. £1,2 ~ Qx/2l/ x x,2 (see below), then we take 
into account the longitudinal momenta in the numerators of effective vertices (I3.3p . 



91 



k 2 



91 k 2 



91 - k'2 
> 

q 2 k\ 
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FIG. 5: Effective vertex in QMRK approach [201 ] . Last diagram with effective 3-gluon vertex drops 
out in projection to the color singlet final state. 

The SU(3) Clebsch-Gordan coefficient (3i,3fc|l) = 5 ik /^W c in Eq. (Q projects out the 
color quantum numbers of the qq pair onto the color singlet state. Factor 1 / y/N~ c provides 
the averaging of the matrix element squared over intermediate color states of quarks. 

Therefore, we have the following amplitude 

(3.4) 

This amplitude can be simplified by using Dirac equations for quark/antiquark spinors 

u\ q {ki)ki = m q ux q {ki), k 2 v Xq -(k 2 ) = -m q v Xq -(k 2 ), k\ = k\ = m 2 q . (3.5) 

Moving k\ to the left, and k 2 to the right, until they disappear upon acting on the spinor, 
we finally get 



2VK d Ux « {kl \ql-2(k l q 1 ) 1 1 ql-2(k 2 q l ) ]V ^ k2} ■ 



(3.6) 



Amplitude of fusion of two off-shell (reggeized) gluons g*g* — > qq turns out to be explicitly 
gauge invariant. Indeed, by direct calculation we see that the gauge invariance over the first 
gluon line is satisfied: 

?TC„ = 0. (3-7) 

Now due to momentum conservation q\ + q 2 = k\ + k 2 we may rewrite the amplitude ( 13.41) 
as follows 

(3.8) 
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Thus, the gauge invariance over the second gluon line is also satisfied: 

Q^xX,^ = °- ( 3 - 9 ) 
Comparing Eqs. (I3.6P and (13 .8p we see that the amplitude is symmetric w.r.t. interchanges 
qi 4-> q 2 and ji -H- v as it should be. 

Taking into account definition (13. 2 p and momentum conservation (12.11) and using the 
gauge invariance properties (13. 7p and (13. 9p . we get the following projection to the light cone 
vectors (so called "Gribov's trick") 

= nln v V^_ „ v = t fi-^-ti-^ v™ = t*L.&Lv** ( 3 . 10 ) 

Last expression shows that an important consequence of the gauge invariance is the vanishing 
of the matrix element of the effective qqqq- vertex between on-mass-shell quark and antiquark 



states in the limit of small gu_ and q 2 ± [17], [22 



V%%->0 for q 1± orq 2± ^0, (3.11) 



The normalization of polarization vectors coincides with that of Ref. [23]. Now using 
Eqs. (I2.3p . (13.1 Op and (13 .3p we finally get the following qq production vertex 

v c lC2 _ 2d 2 , C1C2 - /, \ ( QuA. - 2(k 1± q 1± ) qiqi±-2(k 2± q 1± ) \ 
v \aXa- ~ ~77^ 7rr 6 u x q [Ki) — -2 , - , — q 2 ± - q 2 ± , n r — v Xq {k 2 ). 



^ M* q± VK ^ V1; V <£l-2(*i?i) ^ ^ q 2 1± -2(k 2gi ) 

(3.12) 

It is interesting to note that this vertex function would be equal to zero if one substitutes 
Qi i-e. when one neglects the longitudinal components of gluon momenta qi/ 2 ,i 

putting x\ — > or x 2 — > 0. So, it turns out that the longitudinal momenta play a critical 
role in diffractive production of qq pair and cannot be neglected. At the same time, we keep 
the gluon virtualities in the propagators in Eq. (I3.12p as they apparently become important 
in the small quark masses and quark transverse momenta. 

It is worth to notice that the mass terms disappear when applying the Dirac equations 
(13. 5p . The quark mass m q is present in the spinors and in the scalar products only. So, we 
see that for massless quarks the production amplitude (I3.12p has the same covariant form 
as for massive ones. 

In both particular cases of p' 1± = —p 2 ±_ and in the forward limit \p' 1± \ = |p' 2 ^| — > 0, we 
have qi± = —q 2 ± = q± and, hence, ki± = —k 2 ± = k±. High-fcj_ jets limit corresponds to 
m q <C k± and q± <C k±. Invariant mass of the qq pair is then given by M^= ~ 4|kj_| 2 , and 
the calculation of the matrix element squared |V^(gi, q 2 )\ 2 for the considered hard subprocess 
(I3.12p in this limit leads to 



where is the relative angle between kj_ and q_i_ vectors. We see now that the amplitude 
(I3.12p in the high-fcj_ limit is not exactly zero, but rather suppressed by a factor ~ q\jk\. 
However, relation (13.131) cannot be used for prediction of the corresponding high-fc^ asymp- 
totics of the diffractive amplitude (13. ip since it contains the hard subprocess amplitude 
V(qi,q 2 ) in the first power integrated over q±. For this purpose we have to consider initial 
expression (13.121) for particular quark helicity configurations separately paying attention not 
only at their asymptotical behavior, but also at symmetry w.r.t. q_i_ <H- — q_i_. Large-/cj_ 
behavior of the 4-particle phase space is important as well. 
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B. QQ center of mass helicity amplitudes for the hard subprocess g*g* — > QQ 

Let us consider now separate quark/ ant iquark helicity contributions of the off-shell gluon 
fusion (hard) subprocess g* g* —> Q\ Q\ q , given by the matrix element in the general covari- 
ant form (153211 . 

The most convenient way is to determine the quark/antiquark helicities in the c.m.s. 
frame of the QQ pair with z axis along the proton beam, so \l\ = — k 2 and k\ 2 = M qq /2. 
For simplicity, we work in the limit of forward scattering, so p' l± = p' 2 ± = 0, so qu_ = 
— Q2± = qo±- In this frame momenta of protons and final-state quarks are 



p?=^(l,0, 0, 1), 



p£=^|(l,0, 0, -1) 



ki = E q (l, 7 sin 6 q cos x, 7 sin 9 q sin x, 7 cos 7 = 
= E g (l, —7 sin 6^ cos x, —7 sin 6 q sin x, — 7cos# g ), 



E„ 



< 1, 



so that the proton and quark energies E li2 , E q and the polar angle of a (anti)quark jet 9 q 
w.r.t. the z-axis are defined as 

P — E\ / q 1 g\ _ ^2 / g gx 

-^g — 2 — */2 Xl ~~ \/2 ' 



7ii+Xj 7 x:j + x^ 



where xf 2 are the Sudakov fractions defined in Eq. (12.101) . The gluon and quark transverse 
momenta (with respect to the proton beam) in the polar coordinates are then defined as 



qo± = g±(cos?/>,sm?/>) 
respectively, and 

\/ 7 2 (x^ + xf ) 2 - {x\ - x\) 



>-i± 



fc^(cos x, sin x), 



fcl 



9 1 q 
x\ + x^ 



2 ry>Q rr>Q 

Kz - ^q q . q- 

l ~ ' l 



kl + k 2 z = E ql .{3.U) 



Using these notations, the different helicity amplitudes V\ \ can be written as follows: 



V, = C 



2 1 k I g j_ ^|k| cos(V> — x) — ik z sin(V> — x) j + M qq k±_ (k z cos(2-^ — 2x) — 



z|k| sin(2'0 — 2x) 



/ 



M qq (k 2 ± + q{ + m\) + 4M qq k ± q ± k z cos(^ - x) - 



2k{q 2 ± (l + cos(2^ - 2x)) + q]_ 



(3.15) 



V + 



-2C e 



,q±m q 



k\ cos(2^ - 2x) + |k| 



AM qq k±q±k z cos(^ - x) - 2fcj_g 2 (1 + cos(2-0 - 2x)) + q]_ 



(3.16) 
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with normalisation C = 2g 2 5 ClC ' 2 / y/N~ c . Up to a phase factor, the helicity amplitudes are 
dependent on the difference ip — x (coming from the scalar product (kj_qoj_)) and thus 
explicitly invariant with respect to rotations of kj_ and qoj_ in the transverse plane or shifts 
of the angles ■0, x — > i/j + 8, x + 5. 

Covariant relations (I3.14p allow us to turn to any desirable frame of reference, in particu- 
lar, to the overall c.m.s. frame, where the (anti)quark longitudinal momentum fractions x\\ 
are define through their rapidities y±, y<i in Eq. (I2.10p and the qq invariant mass is given by 

M qq = ^2mi(l + cosh( Z/l - 2/2 )). (3.17) 

Let us now investigate the hard subprocess g*g* — > qq matrix elements (13.151) and (13. 16[) in 
two limits of high and low-fcj_ jets separately. 

Physically interesting case is for qq dijets with very high invariant mass M qq ^> m q , 
where the KKMR QCD mechanism [4| based on the /c^-factorization and Sudakov evolution 
is strictly justified. If one looks at centrally produced jets yi ;2 —> 0, then according to 
Eq. (13.1 7p the only way to produce the large invariant mass M qq is to consider high-fcj_ jets 
limit \kx\ 3> m q , \q±_\. As follows from Eq. (13.141) . such a limit corresponds to the quark and 
antiquark at central rapidities yi^ ~ 0, x\ ~ x\ and the qq invariant mass M qq ~ 2k±. 

We see from Eqs. (13.161) that V ++ — > in the quark massless limit m q — > as it should 
be, so it is generally suppressed with respect to V+_ in the limit of large invariant mass M qq 
and high-fcx quarks, i.e. for M qq , kj_ 3> m q , q±. Indeed, in high-fc^ quarks limit, the helicity 
amplitudes are 

2 2 

V+- ~ -iC^-sin(20-2x), V ++ ~ -C e^^ 1 (1 + cos(20 - 2x)). (3.18) 

Then summing up the non-zeroth contributions |VA g A ? | 2 we easily recover Eq. (I3.13P for 
unpolarized hard matrix element squared in the considered limit. However, V + - actually 
drops out in the integration over ijj in the diffractive amplitude due to antisymmetry with 
respect to the interchange ip <H- —ip. So the diffractive cross section for the qq pair production 
is given by V ++ contribution only and gets significantly suppressed in the high-Zc^ limit by 
a factor ~ q\m q jk\. So, in this case quark masses m q and off-forward corrections become 
important. This is in agreement with the common belief that the Higgs CEP background 
is supposed to be small in very forward and quark massless limits for centrally produced bb 
jets (with small rapidities), and agrees well with the J z = selection rule (lij . 

However, the particular high-fc^ limit does not explain the whole story. Observation we 
have made above means only that quark high-Zc^ contributions in the central rapidity region 
may be strongly suppressed with respect to \ow-k± jets, where the gluon transverse momenta 
q±_ and quark masses are significant, but it does not tell us that the whole bb background 
for Higgs production is suppressed. 

In our previous analysis of the exclusive open charm production in Ref. [l2| it was shown 
that the dominant contribution to the cc dijet cross section comes from relatively small 
quark transverse momenta \k±\ ~ 1 GeV. The same should hold for bb CEP relevant for 
the Higgs background. Indeed, resolving relation (I3.17P with respect to typical rapidity 
difference \yi — y<i\ = Ay, neglecting quark transverse momenta \k±\ m q and keeping 
only the 6-mass contributions rrib — 4.5 GeV at fixed M qq ~ Mh = 120 GeV, we get 
Ay ~ 6.6, x\ ^> x\. So the low-fcj_ quark jets can provide a contribution to irreducible 
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background for Higgs CEP, if their rapidities are yi : 2 — ±3.3. Then the invariant mass of 
the quark/ ant iquark pair is given by their longitudinal fractions (or rapidities) only 

M qq ~ 2m gV / cosh(?/ 1 )cosh(?/2), \yi,2\ > 1 • (3.19) 

In the last kinematical situation the helicity amplitudes ( I3.15p . ( I3.16P are not suppressed by 
a large denominator, and significant contributions can be obtained. Indeed, 



y + _ „ ~" 9 -L y ++ ^ Le ^ (3 20) 

2 v/coshQ/i) cosh(y 2 ) m q {m 2 q + q\) ' ++ 2^cosh(|/i) cosh(y 2 ) m l + Q± 

which in the considered limit q± ~ k± -C m q up to a common phase factor coincide with 



our previously published result in Ref. [11] . Again, analogously to the previous case V + - 
drops out in the integration over if) in the diffractive amplitude. The leading symmetric 
w.r.t. ip -H- — ip contribution to V-\ 

M qq m\ 

is again extremely suppressed w.r.t. V++ in the considered low-fc^ asymptotics. 

Numerically, dominant low-fc^ contribution coming from V ++ amplitude (of course, at 
not extremely large 2/1,2) in the case of 6-jets can lead to a dominant contribution to the 
exclusive background for the Higgs CEP. In this asymptotics, the quark mass m q plays an 
important role since it comes into the denominator in Eq. (13. 20 p . Precise evaluation of the 
corresponding signal, however, demands employing the formulae for the hard amplitudes 
in the general form given by Eqs. ( 13. 15ft and ( I3.16p . Detailed numerical investigation of 
contributions from different parts of the phase space will be presented below in the Results 
section. 



C. Unintegrated gluon distributions 

In our approach we use unintegrated gluon distributions as proposed by Khoze, Martin 
and Ryskin (see e.g. Q). In Ref. 0] slightly different unintegrated gluon distributions taken 
from the analysis of Ivanov and Nikolaev [24[ were used. These gluon distributions have 
been adjusted to the deep-inelastic HERA data. In addition the authors have shown that 
their off-diagonal UGDFs provide a good description of the Tevatron data on exclusive 
dijet production [§]. The KMR UGDFs discussed in the present paper also reasonably well 
describe the dijet data [Io|, 25] . 



Let us now consider in detail the couplings of gluons to protons. At the parton level, 
we assume that hard active gluons (carrying the momentum fractions £1,2) an d screening 
gluons (carrying the momentum fraction x' ^1,2) couple to a quark line in the pro- 
ton in the normal way. In order to turn to the hadron level, the factor CVa s (/ig oft )/7r 
23[ is absorbed into the off-diagonal unintegrated gluon distribution function (UGDF) 



fg S (x', xip, qf/2±i 5oJ-> ^f! t)- The absorbed coupling a s (fi 2 o{t ) corresponds to the coupling of 
the screening gluon with virtuality /ig 0ft ~ <2o ± to a quark in the proton, whereas the coupling 
of the active gluons to the qq central system or to a quark in the proton is purely perturbative 
(given at the hard scale fxp ~ M qq -±) and enters to the hard subprocess amplitude. 
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In the forward limit the following factorization is assumed 



fg (%', xi )2 , q 1/2 ±, ? ±> Vf'i *) = fg ( x '> x i,2> 9i/2±> ?o±> exp(6t/2) , (3.21) 

with the slope parameter b ~ 4GeV -2 j26|. In the x' <C xi^ limit the off-diagonal UGDFs 
can be written as 0, H 



fg 0', zi,2, g 1/2± , g ±. M - R g( x> ) ■ fg( x i,2, ?i/2j.» Mf)> (3-22) 

where the skewedness parameter i? 9 ~ 1.2 — 1.3 is roughly constant at LHC energies, which 
accounts for the single logQ 2 skewed effect [2^] and gives only a small contribution to an 
overall normalization uncertainty. 

Another more symmetrical prescription for skewed UGDFs was introduced in Refs. [3(J 



3 if ] . It is inspired by the positivity constraints for the collinear Generalized Parton Distri- 
butions j32| , and can be considered as a saturation of the Cauchy-Schwarz inequality for the 
density matrix [33]. It allows us to incorporate the actual dependence of the off-diagonal 
UGDFs on longitudinal momentum fraction of the soft screening gluon and its transverse 
momentum in explicitly symmetric way: 

fi/2 g - \] fg( x i,2, q 2 1/2± , /4) ■ fgfa', 9^ /4>ft) > x ' ~ ^ • ( 3 - 23 ) 

As we see it explicitly depends on x' ~ Qq±/s. It works well in the description of the recent 
CDF data on the central exclusive charmonia production (3lL and the precise HERA 
data on the diffractive structure function [35] . Model (13. 23j) implies the factorization of the 
generalized UGDF into the hard part depending on a hard scale \xf and X\p describing the 
hard gluon coupling to the proton, and the soft part defined at some soft scale // so ft and 
small x' X\ t 2- Together with the factorization in transverse momentum space, model 
( I3.23P provides the QCD factorization of the diffractive amplitude in the full momentum 
space. 




FIG. 6: Gluon densities as a function of longitudinal momentum fraction x at the scales Q 2 = 1, 2 
and 5 GeV^given by the global parameterizations CTEQ6L1 [44], GRV94HO [42], MSTW2008LO 
and NLO 0. 



In the considered kinematics the diagonal unintegrated densities can be written in terms 
of the conventional (integrated) densities xg(x, qj_) as [28 



f g (x,ql,fi 2 



d 

d In q\ 



[xg(x,qi)JTg(ql,fj, 



(3.24) 



12 



where T g is the Sudakov form factor which suppresses real emissions from the active gluon 
during the evolution, so that the rapidity gaps survive. It is given by 



T g (q. 



exp 



^ dk 2 t 



a s {k] 



2n 



l-A 



dz 



(3.25) 



where the upper limit is taken to be 



A 



k 



k± + aM t 



(3.26) 



The KMR group used a = 0.62 36]. It was argued recently that a = 1 should be used 
instead |37|], so in numerical calculations below we adopt the last choice (for these two 
choices the final results for the Higgs CEP cross section differ by about a factor of 2, which 
is not negligible). In addition, following Ref. |37| we use the factorization scale fip — M, 



qq 

as compared to the KKMR choice /z™ R = M q q/2. We will discuss the sensitivity to the 
factorization scale choice below when presenting numerical results. 

In general, employing diagonal UGDF in the form f)3.24p we encounter a problem of 
poorly known gluon PDFs at rather low 2 12 and especially small gluon virtualities qj_. For 
an illustration of the corresponding uncertainties, in Fig. [6] we show several gluon PDFs as 
functions of fraction x at evolution scale ~ q\ fixed at values 1, 2 and 5 GeV 2 characteristic 
for the exclusive production of Higgs boson. We see that at x < 10 -3 the PDF uncertainties 
may strongly affect predictions for not sufficiently large gluon transverse momenta. In this 
sense, the precise data on the diffractive and central exclusive production could be used for 
making constraints on the PDF parameterizations (35| . 

Testing other models of UGDFs, different from Eq. (13.241) . may be important for estima- 
tion of an overall theoretical uncertainty of our predictions and their stability. 



IV. ELECTROMAGNETIC 7* 7 *-FUSION PROCESS 

It is instructive to estimate the QED contribution to the central exclusive bb production 
illustrated in Fig. [3 



Pi ^ 




FIG. 7: The QED 7*7* fusion mechanism of the exclusive qq production. 
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In the forward limit of small momentum transfers ti t2 ( \U 2 1 -C 4m^) the matrix element 
for pp —y ppqq reaction via 7*7*-fusion can be written as [30(] 

M 1 * 1 * * eF x {t x ) iPl+ /' r V^( gi , q 2 ) ^ + /* r efifo) , (4.1) 
t\ t 2 

where Fj {U ) and Fi(t 2 ) are Dirac proton electromagnetic form factors, and the 7*7* — > qq 
vertex [31| has analogous form as (13. 4p . i.e. 

(4.2) 

Momentum conservation dictates us the following decompositions of the photon momenta 



into the longitudinal and transverse parts w.r.t. c.m.s. direction [31 



qi = X1P1 + f -jP2 + qi±, Q2 = X2P2 + jPi + ?2±, 9?/ 2 ± - *i,2(l - £1,2) , (4.3) 
where t x>2 = ?i 2 - Due to gauge invariance we have (similarly to (I3.10P ) 

V^*(quq 2 ) = (pi+tiTV^iquq^+tir = ^^(ft,®,). (4.4) 

so the matrix element squared |^/ 7 * 7 * | 2 is proportional to the gluonic one found in Eq. (13. 13ft . 
The photon virtualities in the relevant limit disappear t\ j2 — > 0, so q± — > leading to van- 
ishing of the 7*7* — > qq amplitude in the high-Zc^ limit and massless quarks (see Eq. (I3.13P ). 
Thus, to estimate 7*7* contribution to exclusive production of quark jets we have to take 
into account subleading corrections in the m g //c^-expansion. This means that the elec- 
tromagnetic mechanism may give some contribution for moderate and large quark masses, 
which will be evaluated numerically in the Results section. 



V. OFF-SHELL EFFECTS IN CENTRAL EXCLUSIVE HIGGS PRODUCTION 

As one can see from Eq. (I3.10p . the subprocess vertex V\ q \ q of g*g* — > bb is proportional 
to gluon transverse momenta squared due to projection ~ <li±<l 2 ±- When considering the 
irreducible bb background for Higgs boson production in the amplitude V\ q \- we can neglect, 
in principle, the gluon virtualities in comparison with quark transverse momenta k± and 
6-quark mass m q as it was done in e.g. Ref. jl]. 

However, as demanded by fc^-factorization framework, it may be instructive to analyze in 
which region of the phase space the gluon off-shell effects may play a significant role (if any), 
and whether it is possible to see such effects in experiment. A complete calculation of the 
off-shell effects in inclusive Higgs boson production was performed in Ref. j38[ (for the Higgs 
boson production in /c^-factorisation, see also Ref. (39|). It was shown there that the off- 
shell effects can significantly affect the distribution of Higgs boson cross section in azimuthal 
angle between fusing gluons in a very close vicinity of = tt/2. The calculations of the 
central exclusive Higgs production rates in the on-shell gluon approximation are well-known 
Let us now investigate off-shell effects in Higgs CEP. 
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A. Off-shell effects in the hard vertex g*g* — > H 



Tensor decomposition of the hard subprocess amplitude g*g* — > H can be written in the 
following general form [38] 



T? u (qi, q 2 ) = iS*~ ([(fcfc)^ - ft^Jd + 



+ 



q\ 



(Q1Q2) 



q2, fi q2, l 



q\ q\q\ 
qi ttl qi,u + 



(qm, 



G 2 ), (5.1) 



where a s is the strong coupling constant, v = (Gpy/2)^ 1 ^ 2 is the electroweak parameter of 
the Standard Model, a, b are the colour indices of two virtual gluons with momenta q±, q 2 . 
Let us introduce the dimensionless parameters 



M 



X 



H 



4m 2 



> 0, 



Xi 



Am 2 



< 0, 



X2 



< 0, 



so the heavy quark limit corresponds to x, Xi> X2 — > 0. In the case of heavy Higgs production 
we have Mfj 3> |qi_i_| 2 , | q 2 _L 1 2 ? so in the expansion of the form factors we have to take into 



account powers of x higher than powers of gluon virtualities xu X2 [3 

Gi(x,Xi,X2) = 



1 7 2 2 11/ 
1 + 30 X+ 21 X + 30 (Xl + X2) + - 



1 4 

G 2 {x,Xi,X2) = - — (x-Xi-X2)-—X 2 + --- 
45 315 



(5.2) 



These expansions will be sufficient for our present calculations. 

Let us turn now to the discussion of the exclusive diffractive Higgs production within the 
KKMR double diffractive mechanism Q. The hard subprocess vertex entering the diffractive 
amplitude (13. ip V = V(qf ± q 2 ± , Pj_) with explicit taking into account gluon virtualities reads 



V^ H (qf ± ql ± , Pi) = n+rCZ^fa, q 2 ) = 1 ^^Z^fa, q 2 ), tfr£ = £1$ = 0, (5.3) 



S Xi x 2 



where the amplitude of the gluon fusion T^(qi,q 2 ) is defined in Eq. (15. ip . Contracting 
indices and introducing the transverse Higgs mass, 



sx lX2 = M 2 H -Pl = M* ± , Pi = 
we get in terms of form factors (15. 2p 

,ab a s |qi±||q2±| 



+ qL + 2|qi±||q 2 ±| cos0), (5.4) 



yab 



45° 



71 



cos (j) G\ 



2M^ Jq 1± ||q 2± | 
(M 2 H + q? ± + q^) 2 



C 5 



(5.5) 



In the limit of real gluons for not extremely heavy Higgs we have asymptotically 



V. 



ab 

gg^H 



-iS*—- qi±q 2 ± • -, 

TT V 3 



lqi±q2±, 



|qi±||q2±| cos0. 



(5.6) 
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Substituting this into the amplitude (13. ip . we get after summation over colour indices 

A A on — shell 1Xa s S / 2 ( ( ll± ( l2±j j. B I I 2 2 r\ eofi ( 12 2 r \ 

M excl dc io± o — 2 o f g ,i{xi,x ,q 0± ,q 1± ,t)f gt2 {x 2 ,x ,q 0± ,q 2± ,t 2 ), 

v o J q 0± q 1± q 2± 

Next-to-leading order contribution in the cross section of the hard subprocess a(gg — > H) 
can be accounted for by a factor K^lo — 1-5 assuming that the NLO corrections factor -K'nlo 
in the q*q* — > H vertex is the same as in the H — >■ gg width provided that \ V gg ^H | 2 ~ T(H — > 
gg) [40]. As for the irreducible Higgs background, the one- loop corrections to gg — > bb were 
calculated in Ref. 0. For simplicity, in the current analysis we do not take into account the 
NLO effects and concentrate only on the leading order contributions in the hard subprocess 
part. 

At high energies the cross section for the exclusive Higgs boson production can be ex- 
pressed as: 

da pp ^ P Hp = | M | 2 • d 3 PS, d 3 PS = -^—dhdhdynd® . 



B. bb signal of the Higgs decay 



In order to estimate the observable signal from the central exclusive production of Higgs 
in the bb decay channel, one can multiply the matrix element squared \M. exc i\ 2 by the rela- 
tivists Breit-Wigner distribution over the invariant mass p = M bb of bb pair with a proper 
normalization 



PbM 



i 



(5.7) 



vr \p? - M 2 H } 2 + [pT™{p)] 2 ' 

and then integrate it out over the 4-particle p(bb)p invariant phas e sp ace. In Eq. (15 .7p the 
total Higgs decay width into fermions and gluons is (see e.g. Ref. [411 ] ) 

r f J(M H ) + t™{m h ) , 



r f J(M H ) 



H 



H 



32ttM^ 



3[m 2 b (M H )+m 2 c (M H )} 



1 + 5.67^^ + 42.74 



7T 



« S (M 



11, 



71 



(5.8) 



+ 



r »(M„) = / M K a - (M ">h + 17.91667 a - (M " ) 



288 ttM 2 , 



71- 



71 



and the partial Higgs decay width into bb channel is given by 



3g 2 M H , 



32 71 M 2 , 



m 2 b {M t 



a s (M H ) 
5.67 V ; +42.74 



7T 



ot B {M i 



H, 



7i 



(5.9) 



Above, g 2 = 0.42502 is the electromagnetic coupling constant, and m c ,b{p) and a s (p) are the 
running quark masses and the QCD coupling constant, respectively. For example, at Mh = 
100 GeV, they are m c (M H ) = 0.542 GeV, m b {M H ) = 2.676 GeV and a s {M H ) = 0.12121. 
Substituting these values into Eq. §EM we S et the total deca y width r^(M fl ) = 2.268 MeV 
which coincides with the number given in Ref. (4lj . 
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VI. NUMERICAL RESULTS 



In this section we will present differential distributions for central exclusive heavy quark 
dijets and Higgs production. We show also distributions of b and b quarks from the decay of 
the Higgs boson. We start from the presentation of the results for heavy quark (cc and bb) 
production. Very important part of the analysis below concerns the bb background below 
the bb Higgs signal. 



A. Differential distributions for exclusive cc and bb pair production 

We start our presentation from the distribution in heavy quark invariant mass. In Fig. [H] 
we show distributions for cc (left panel) and bb (right panel) for different gluon collinear 
distributions used to generate UGDFs, for quark rapidities y q ,y q G (—2.5,2.5). We show 
distributions for QCD diffractive mechanism as well as for the QED 7*7* fusion. Only at 
small invariant masses is the distribution for cc higher than that for bb. The position of the 
peak of the distributions depends on the quark mass and is placed slightly above 2m q , both 
for diffractive and QED components. Relative contribution of QED mechanism grows with 
increasing invariant qq mass, and starts to dominate at M qq > 60 GeV for cc and M qq > 120 
GeV for bb. 




M cc (GeV) M (GeV) 



FIG. 8: Invariant mass M qq distributions for the central exclusive production of cc (left panel) and 
bb (right panel) dijets at LHC. Results for diffractive component are shown for a few GDFs from the 
literature: GRV94 NLO (dashed line), GJR08 NLO Q (dotted line), CTEQ6 NLO (solid 
line) and MSTW2008 NLO (long-dashed line). QED 7*7* fusion contribution is illustrated by 
the dash-dotted line. Experimental cuts on quark (antiquark) rapidities \y q \ < 2.5 are included. 
The distribution is then integrated over whole invariant qq mass range 2m q < M qq < 200 GeV in 
order to get the total cross section. 



There is a strong sensitivity of the invariant mass distribution on the gluon PDFs. In 
particular, the difference of the results for the GRV94 NLO [II] and CTEQ6 NLO 0] is up 
to an order of magnitude in the peak. The QED contribution is found to be smaller than 
that for the diffractive mechanism but it is very important as background to the exclusive 
Higgs boson production. In this calculation only the Dirac F\ proton electromagnetic form 
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factor is included. The contribution of the Pauli F2 proton electromagnetic form factor is 
expected to be negligible. 

For the bulk of the cc or bb production M q q is only slightly larger than the quark masses. 
Nevertheless what matters is that M qq - is in the perturbative region. We can go down with 
transverse momenta of gluons as low as q\ ~ 0.4 GeV 2 . Of course it is not easy to predict 
what happens in the nonperturbative region. In the present paper we concentrate on large 
quark- ant iquark invariant masses where the issue is not so important. 




FIG. 9: Invariant mass distributions of centrally produced cc pair (left panel) and bb pair (right 
panel) for different (anti)quark helicities \ q \q = ++, — V and for the sum over all quark helicity 
states. CTEQ6 PDF was used here. Kinematical constraints are the same as in Fig. [8l 

There was recently [7] a discussion about the contribution of different helicity states to 
the cross section. In Fig. [9] we show individual contributions for different quark helicities 
W = ++ and — h (other helicity contributions are the same due to symmetry as discussed 
above). The contributions of the same and opposite quark helicities are rather similar in 
the broad range of the quark-antiquark pair invariant masses. In particular, they are almost 
identical in the region of typical light Higgs mass. 

The cross section for the diquark production strongly depends on the quark masses thus 
underlying the importance on the finite 6-quark masses for the CEP Higgs background 
evaluation. This is encoded in the matrix elements discussed in the theory section. In 
Fig. [10] we have collected the results for the total (left panel) and differential in quark- 
antiquark invariant mass (right panel) cross sections. The total cross section slightly grows 
with the quark mass. The growth is, however, slower than in Ref. 0] where the matrix 
element is proportional to the quark mass. Taken a typical misidentification probability and 
the fact that light quark cross sections are smaller than that for bb in the Higgs region, the 
latter is in practice the only troublesome background. Our result is very interesting in the 
context of diffractive dijet production. Recently the CDF collaboration has measured the 
corresponding cross section at the Tevatron {3J. The quark jets are usually neglected and 
only gluonic jets are included in theoretical calculations (for the gluonic jets analysis, see 
e.g. Ref. [9J). The quark jets contribution to the CDF data will be discussed elsewhere (2o| . 

Let us come now to the rapidity distributions. The distribution in quark (antiquark) 
rapidity in the detector interval, integrated over whole invariant mass range of M qq - 6 
(2m q , 200) GeV, is shown in Fig. [TTJ The distribution for charm quarks is flatter than 
that for the bottom quarks. For comparison, we also show the corresponding distributions 
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FIG. 10: The EDD qq cross section in the vicinity of the Higgs mass as a function of quark 
mass (left panel) and invariant mass distributions of the EDD qq production for different quark 
masses. Kinematical constraints are the same as in Fig. [51 



for different combinations of quark helicities \ q \q = ++ and — h As for the quark-antiquark 
invariant mass distributions, the same and opposite helicity contributions are similar. The 
rapidity distributions shown are dominated by the low quark-antiquark invariant masses. 
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FIG. 11: Differential distributions in quark (antiquark) rapidity of EDD c or c (left panel) and b 
or b (right panel) . Kinematical constraints are the same as in Fig. [8j 

In Fig. [12] we show distribution in the rapidity of the qq pair. Here the distribution of the 
opposite quark helicity is much natter than the distribution of the same quark helicity. This 
distribution may, however, be slightly biased by the limitation of the individual rapidities 
of the quark and antiquark. 

Let us come now to transverse momentum distributions. In Fig. [T3] we show distribu- 
tion in quark (antiquark) transverse momenta. These distributions are extended to large 
transverse momenta with the peak at about 2 GeV for c/c and 3 GeV for b/b. This is fully 
perturbative effect and is encoded in the g*g* — > qq matrix elements discussed in Sections 
III and IV. One can clearly see the dominance of the opposite sign helicities contribution at 
large transverse momenta. For comparison, we show also distribution in proton transverse 
momenta (dashed line). In contrast to quarks (antiquarks) pj_-distributions, they are order 
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FIG. 12: Differential distributions in rapidity of the cc pair Y cS (left panel) and bb pair Y h i (right 
panel) . Kinematical constraints are the same as in Fig. El 




FIG. 13: Differential distributions in transverse momenta of quark (antiquark) (solid line) and 
outgoing protons (dashed line) of EDD cc (left panel) and bb (right panel) cross section. Kinematical 
constraints are the same as in Fig. [HJ 



of magnitude narrower and concentrated below 1 GeV with maximum at about 0.3 GeV. 
These distributions are controlled by a nonperturbative proton form factor and are thus 
sensitive to internal structure of the proton. 

The distribution in the total transverse momentum of the qq pair |p g g±| (by definition 
Pgg± — Pq± + Pq±) is shown in Fig. [TH It is much narrower than that for the individual 
quark (antiquark). The maximum of the cross section is at about 0.5 GeV. Similarly as the 
distributions in the proton transverse momentum this distribution is fully nonperturbative 
and related to the slope of the nucleon form factors. 

Finally, let us turn to azimuthal angle correlations. We will consider correlations between 
outgoing quark jets, as well as between outgoing protons. 

In Fig. [15] we show correlations between outgoing jets without extra cuts on jets trans- 
verse momenta. Even without such cuts the quark and antiquark are strongly correlated 
with a preference for the back-to-back configuration. The deviation from the back-to-back 
configuration is caused by the transverse momenta of gluons in the ladder. If there were no 
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FIG. 14: Differential distributions in transverse momenta of quark-antiquark pair |p g gj_|. Kine- 
matical constraints are the same as in Fig. [8l 




FIG. 15: Differential distributions in azimuthal angle between c, c jets (left panel) and b, b jets 
(right panel) . Kinematical constraints are the same as in Fig. [HJ 



transverse momenta of initial gluons, final jets would be back-to-back which follows from 
the kinematics of the process. There is a stronger helicity correlation for the opposite quark 
helicities than that for the same quark helicities. The correlation would even increase when 
imposing extra cuts on quark (antiquark) transverse momenta. 

In Fig. [16] we show correlations between outgoing protons. In contrast to quarks (anti- 
quarks) protons are almost decorrelated. This can be understood post factum taken com- 
plicated gluonic ladders spanned between protons and quarks. The soft rescattering effects 
could further modify the distribution (see e.g. 26, 46|). 



Finally, we would like to show a two-dimensional distribution which is very useful when 
discussing background to the exclusive Higgs boson production in the bb channel. In Fig. [T7] 
we show the distribution in M b i and transverse momentum of the quark (pb±) for EDD 
(left panel) and QED (right panel) mechanisms. One can clearly see that a fixed mass 
(e.g. mass of the Higgs) can be obtained both for high and low transverse momenta of the 
quark jets. The latter case can be realized when the quark rapidities are large. For EDD 
contribution one can remove such cases by imposing cuts on jet transverse momenta. One 
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FIG. 16: Differential distributions in angle between outgoing protons corresponding to cc EDD 
(left panel) and bb EDD (right panel) . Kinematical constraints are the same as in Fig. [8l 



EDD pp^ppbb CEP yy^bb 




FIG. 17: Two-dimensional distribution in invariant bb mass M b i and (anti)quark transverse mo- 
mentum pb± . Kinematical constraints are the same as in Fig. El 



could equivalently limit quark (antiquark) rapidities. We will return to these correlations 
when discussing the Higgs background. 



B. Central exclusive production of Higgs boson 

Due to a very large hard scale of the process ~ M^, the influence of typically small 
gluon virtualities in the amplitude of hard subprocess amplitude (15.51) . as well as the role of 
form factor G 2 , in the exclusive diffractive Higgs production turned out to be quite small, 
in analogy to the inclusive case [38| . 

The integrated cross section of diffractive Higgs production at the LHC energy = 14 
TeV, taking into account the "effective" gap survival factor (S 2 ) ~ 0.03 [4], calculated for 
typical Higgs mass Mh = 120 GeV is o to t ^ 1 fb. Our result is smaller than that found by 
the Durham group. As discussed above this is mostly due to different choice of the scale 
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of the Sudakov form factor. The result of Cudell, Hernandez, Ivanov and Dechambre 0] is 
closer to our result but still slightly bigger. This is probably due to different unintegrated 
gluon distribution. In particular, the Ivanov-Nikolaev UGDF used in their analysis includes 
also a nonperturbative piece fitted to the data. 
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FIG. 18: Total cross section for exclusive Higgs production for different gluon PDFs from the 
literature. The calculation was done including off-shell effects. 



In Fig. [18] we show the total cross section for exclusive production of Higgs boson as a 
function of the Higgs mass for different gluon distributions for yfs = 14 TeV. The difference 
between different gluon PDFs comes mainly from a different lower cut-off parameter for gluon 
transverse momenta in different gluon distributions. This is necessary and is dictated by the 
construction of different UGDFs. In particular, different groups choose different initial scale 
for QCD evolution and going below it often leads to unphysical solutions (negative glue for 
instance). This forces one to put lower cut-off at the value of the initial scale. The cross 
section for exclusive Higgs production obtained here is rather small 2 . 

We have made the calculation of the cross section in the limit of real gluons in the hard 
part (15. 6 p (cr|f), as well as with an account of gluon virtualities (15 .5p (cr£f )• Contribution 
of non-zeroth qf, q\ in form factors G\ % i turns out to be negligibly small; difference between 
o°ft and cr^f is formed mainly by the second form factor G 2 , and gives about 6 %, so it is 
much smaller than other theoretical uncertainties of the approach. The overall uncertainty 
of + Higgs CEP cross section was estimated in Ref. |4| to be up to a factor of 2.5. 

In Fig. [in] we show a two-dimensional distribution of the Higgs in its rapidity and trans- 
verse momentum. The Higgs production is concentrated around rapidity y = and the 
cross section quickly drops with Higgs transverse momentum. In Fig. [20] we show respective 
projections on rapidity (left panel) and transverse momentum (right panel). The maximum 
of the transverse momentum dependence occurs at about 0.4 GeV. The distribution reflects 
a convolution of the nucleon form factors, i.e. is of purely nonperturbative nature. 

Finally, we focus on angular correlations (see Fig. I2~lj) . In the figure we show distribution 
in azimuthal angle between outgoing protons. As for the exclusive production of heavy 
quarks there is a very small correlation between outgoing protons. 

Note that the distribution in relative azimuthal angle between protons PP strongly differs 



from the distributions in azimuthal angle 



y qiq2 



between interacting gluons ~ cos 



y qiq2 



due 



Similarly small cross sections have been obtained very recently 15[ when this paper was already finished. 
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EDD pp^pHp 




FIG. 19: A two-dimensional distribution in the Higgs rapidity (y) and Higgs transverse momentum 
(p_i_) for CTEQ6 gluon PDF. 




FIG. 20: Rapidity (left) and transverse momentum (right) distributions of the Higgs boson. 
CTEQ6 PDF was used in this calculation. 



to the convolution with momentum transfer q ± of the screening gluon. Only in the case 
when Higgs boson production is governed by e.g. pomeron-pomeron (or 7*7*) fusion, the 
angle between pomerons (photons) coincides with the angle between outgoing protons, and 
corresponding distribution has cos 2 pp -dependence. This fundamental difference of the two 
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mechanisms was observed also in other processes 

We have not discussed yet the influence of the off-shell effects in the matrix element 
on differential distributions. In the limit of real gluons the differential distributions are 
practically unchanged, so the corresponding off-shell effects are hard to observe taken the 



other theoretical uncertainties. Moreover, unlike for inclusive production case 38], shapes 
of differential distributions of Higgs CEP are not sensitive to the off-shell effects since they 
get averaged out effectively when the off-shell matrix element (15.51) is integrated over q j_ in 
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FIG. 21: Differential distribution in angle between outgoing protons for central exclusive Higgs 
boson production. CTEQ6 PDF was used in this calculation. 

the diffractive amplitude (13. ip . 



C. Irreducible 66 background for exclusive Higgs production 
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FIG. 22: The 66 invariant mass distribution for yfs = 14 TeV and for 6 and 6 jets from Higgs decay 
in the rapidity interval —2.5 < yb < 2.5 corresponding to the ATLAS detector. The absorption 
effects for the Higgs boson and the background were taken into account by multiplying cross section 
by the gap survival factor (S 2 ) = 0.03. The left panel shows purely theoretical predictions, while 
the right panel includes experimental effects due to experimental uncertainty in invariant mass 
measurement. The left peaks (bumps) correspond to the Z° contribution and the right ones to the 
Higgs contribution. 



Now we turn to the analysis of the bb continuum as a background for the bb Higgs signal. 
In the left panel of Fig. [22] we show contributions of several CEP mechanisms to the bb quark 
invariant mass distribution. The diffractive bb and Higgs contributions were calculated for 
a selected (CTEQ6 j44J) collinear gluon distribution. The QED mechanism is also shown 



by the short-dashed line. Natural decay width, calculated as in Ref. [41J, was assumed in 
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this calculation, see the sharp peak at M bb = 120 GeV (assumed arbitrarily for illustration) 
which is not excluded at present by the Higgs searches at LEP [47j and Tevatron (48| . 

As was already mentioned above, the phase space integrated cross section for the Higgs 
production, including absorption effects with (S 2 ) = 0.03 is less than 1 fb which is somewhat 
smaller than that predicted by the Durham group 0]. The main reason is different choice of 
the scale in the Sudakov form factor fj, 2 = M\ instead of [i 2 = as used by the Durham 

group. The first choice was advocated recently by theoretical studies in Ref. 37}. The result 



shown in Fig. [22] includes also the branching fraction BR,(H — > bb) ~ 0.8 and the rapidity 
restrictions. The second much broader Breit-Wigner type peak corresponds to the exclusive 
production of the Z° boson with the cross section calculated as in Ref. [49]. The exclusive 
cross section for -y/i = 14 TeV is 16.61 fb including absorption (28.71 fb without absorption 
effects). The branching fraction BR(Z° — > bb) w 0.15 has been included in addition. In 
contrast to the Higgs case the absorption effects for the Z° production are much smaller 
[49l | . The sharp peak corresponding to the Higgs boson clearly sticks above the background. 
In the above calculations we have assumed an ideal no-error measurement. 

In reality the situation is, however, much worse as both protons and the b and b jets are 
measured with a certain precision which automatically leads to a smearing of experimental 
distribution in M bb . Much better resolution can be obtained by measuring missing mass than 
from the direct measurement of heavy quark (antiquark) jet momenta. In the following in 
spite we will present distribution in M bb (it will mean experimentally distribution in missing 
mass (Mp P )). The two are identical when there are no errors on kinematical variables. While 
the smearing is negligible for the background, it leads to a significant modification of the 
Breit-Wigner peaks, especially of the sharp one for the Higgs boson. In the present paper 
the experimental effects are included in the simplest way by a convolution of the theoretical 
distributions and experimental resolution function 

da exp , . , s f . da 



r d a th 



where the experimental resolution function is taken as the Gaussian function 

G{x) = -^exp(£) , (6-2) 
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with a = 2 GeV, which realistically represents the experimental situation [50j, |51] and is 
determined mainly by the precision of measuring forward protons. In the right panel we 
show the invariant mass distribution when the invariant mass smearing is included. Now the 
bump corresponding to the Higgs boson is below the bb background. With the experimen- 
tal resolution assumed above the identification of the Standard Model Higgs looks rather 



difficult. The situation for some scenarios beyond the Standard Model may be better [52 

Below we wish to discuss how to improve the situation by imposing extra cuts. Before we 
establish how to impose cuts let us consider first a few two-dimensional distributions which 
may help to come to the final solution. 

Let us start from two-dimensional distributions in proton transverse momenta. In Fig. [23] 
we show distributions for the diffractive (left panel), photon-photon (middle panel) and for 
b and b from the Higgs boson decay (right panel) contributions. In the case of the EDD 
and QED continua we are limited to a very restrictive range of invariant masses (117.6 GeV 
< M bb < 122.4 GeV) around the chosen Higgs mass in order to facilitate a comparison of the 
signal and background. While the distributions for the diffractive bb continuum and Higgs 
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FIG. 23: Two-dimensional distributions in outgoing proton momenta pi t 2± of the bb EDD (left) 
and QED (middle) continua integrated in the window 117.6 GeV < M bb - < 122.4 GeV and Higgs 
CEP signal (right panel) . Kinematical constraints are the same as in Fig. |2"21 

are rather similar the distribution for the photon-photon production differs considerably. 
While the first two ones are peaked at sizeable transverse momenta of about 0.3 GeV, the 
photon-photon contribution is peaked at extremely small proton transverse momenta due to 
photon propagators. Cutting off extremely small proton transverse momenta would allow 
to get rid of the photon-photon contribution to a large extent. It is not completely clear 
if this can be done easily experimentally. A Monte Carlo study including the experimental 
apparatus seems to be required. 
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FIG. 24: Two-dimensional distributions in b and b rapidities integrated in the window 117.6 GeV 
< M b i < 122 A GeV for the diffractive QCD background (left panel), 7*7* contribution (middle 
panel) and Higgs CEP (right panel) . Kinematical constraints are the same as in Fig. [22j 

Next let us consider two-dimensional distributions in rapidities of the quark (y b ) and 
antiquark (y b ). In Fig. |2U as in the previous case, we show separately distributions for 
diffractive continuum (left panel), photon-photon continuum (middle panel) and Higgs (right 
panel) contributions. The problem of the background subtraction looks here fairly favorable. 
In the case of the bb continuum production the cross section is maximal when quark and 
antiquark have opposite rapidities at the edges of the main detector (ATLAS, CMS). This 
is completely different for the Higgs contribution where the maximum occurs when y^, yi ~ 
0. Two windows suggesting how to get rid of the major part of the background are shown 
in Fig. [2U the square marked by the dashed line and the area between two parallel lines at 
45°. The consequences of such cuts will be discussed in the following. 
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The situation can be also quantified in a one- dimensional plot in a function of the dif- 
ference of the quark and antiquark rapidities (see Fig. [25]) . The distributions for the signal 
and background are very different. Imposing a cut on ydiff can significantly improve the 
signal-to-background ratio. 




FIG. 25: Distribution in the difference of the quark and antiquark rapidities. Please note an extra 
cut on the bb invariant mass. Kinematical constraints are the same as in Fig. [22] 

In Fig. [26] we show the distribution in the 6-quark rapidity from Higgs decay and from 
a narrow region of bb invariant mass (given in the figure) for the diffractive bb and photon- 
photon components. While the Higgs contribution is concentrated at ?/& ~ the diffractive 
component has maxima at the edges of the central detector. The 77 contribution is rather 
flat across the range of the central detector. The different distributions in the 6-quark 
rapidity of the different components suggest that limiting to midrapidities (i.e. not using 
the whole range of the detector) may help in improving the signal-to-background ratio. 

Further useful handles to improve the situation are the jet transverse momenta which 
can be measured in the central detector. The importance of the cuts on the jet transverse 
momenta is illustrated in Fig. [27] Again we show the three components. While the signal 
(Higgs) contribution is peaked at the transverse momenta being half of the Higgs mass, the 
background contributions are flat or even have local maxima at low transverse momenta. 
Imposing therefore a lower cut on jet transverse momenta can again significantly improve 
the signal-to-background ratio without losing too much of the signal itself. Also, from 
experimental point of view the b (b) jets can be well identified only above a certain cut on 
their transverse momenta. 

Now we wish to quantify the effect of cuts on the bb invariant mass (missing mass exper- 
imentally) distribution. We shall impose cuts in order not to loose too much Higgs signal. 
In Fig. [28] we show the results for several scenarios (cuts). Here we omit the Z° contribution 
and concentrate solely on the Higgs signal. In the left upper corner we show result with 
the cut only on quark and antiquark rapidities (the square in Fig. 121]) i.e. not making use 
of the whole coverage of the main LHC detectors. The signal is now above the diffractive 
background. We also show, by the thin dashed line, the photon-photon background which is 
only slightly smaller than the diffractive one. In the upper right corner we show the result 
for the cut on the quark and antiquark rapidity difference (see parallel thick solid lines in 
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FIG. 26: Distribution in the quark/antiquark rapidity. Please note an extra cut on the bb invariant 
mass. Kinematical constraints are the same as in Fig. [22j 




FIG. 27: Distribution in the jet transverse momentum for different components. Please note an 
extra cut on the bb invariant mass. Kinematical constraints are the same as in Fig. [22j 



Fig. 124)) . The signal-to-background ratio is here similar, except that the cross sections are 
larger. In the lower left corner we show the situation with the lower cut on both quark and 
antiquark jets. The situation is similar as for the rapidity cuts. In order to eliminate the 
photon-photon contribution in the lower right corner we impose in addition a lower cut on 
proton transverse momenta. Now the signal clearly sticks above the background and the 
contribution of the photon-photon continuum is negligible. The cross section for the Higgs 
boson with the cuts is only 2-3 times smaller than that without the cuts. 

The simultaneous inclusion of cuts on quark (antiquark) rapidities and transverse mo- 
menta does not improve further the signal-to-background ratio. It is enough in practice to 
include only one of them depending on experimental convenience. Why it is so is discussed in 
Fig. [291 where we show two-dimensional distributions in 6-quark rapidity and transverse mo- 
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FIG. 28: The bb invariant mass distribution for y/s = 14 TeV and for b and b jets from Higgs 
decay with different limitations in the (anti)quark rapidity y q , (anti)quark transverse momenta and 
outgoing proton momenta p\ 2_i_. 



mentum for the bb invariant mass window around the Higgs mass. Here three well separated 
maxima are seen. One can extract them either by cuts on quark/antiquark rapidities or 
by cuts on quark/antiquark transverse momenta. This explains the equivalence of the cuts. 
Unlike the corresponding distribution not restricted to M b i around Higgs, the distributions 
with the restriction have more complicated structure. While for the Higgs decay 6-quarks 
are produced predominantly at midrapidities with transverse momenta p± ~ Mh/2. The 
two-dimensional distributions for the diffractive continuum and the 7*7* fusion subprocess 
have three maxima: one as for the Higgs decay (smaller) and two other with much smaller 
transverse momenta and larger rapidities (bigger). For diffractive bb continuum the maxima 
at small p± and large rapidities are much bigger. Identification of the maxima experimen- 
tally would be a confirmation of the present predictions. Imposing appropriate cuts in either 
rapidity or transverse momentum would allow to get rid of a large fraction of the diffractive 
background. 

The equivalence of the cuts in rapidity and transverse momentum can be even better 
understood in the two-dimensional distribution in 6-quark transverse momentum and the 
difference between quark and antiquark rapidities. Fig. 1301 shows that for a fixed narrow 
interval of the bb invariant mass the rapidity difference between b and b and the trans- 
verse momentum of b or b are strongly correlated. This is of purely kinematical origin (see 
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FIG. 29: Two-dimensional distribution in quark rapidity and quark transverse momentum for 
different mechanisms. 
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FIG. 30: Two-dimensional distribution in the quark-antiquark rapidity difference and quark trans- 
verse momentum. 
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FIG. 31: The two-dimensional distributions in £1 and £2 for diffractive bb (left panel), 7*7* — > bb 
(middle panel) and CEP of Higgs (right panel) . 

Eq. ( I3.17P ) and demonstrates that imposing cuts on one of the two variables is equivalent 
and completely sufficient. 

Above we have considered proton transverse momentum cuts. Can other outgoing proton 
variables be useful to improve the signal-to-background ratio? In Fig. [31] we show the 
distribution in proton longitudinal momentum fraction losts £1 and £2 defined as: 

£1 = 1 — x^(proton 1 ) , £2 = 1— a;i?(proton 2 ) , (6.3) 
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where z^'s are the Feynman variables of outgoing protons 1 or 2. Only slightly different 
distributions for exclusive Higgs production (right panel), bb EDD (left panel) and QED bb 
(middle panel) continua can be seen. Imposing cuts on £1 > £ cut or £ 2 > £ C ut could slightly 
improve the signal-to-background ratio but at the expense of severe deteriorating the statis- 
tics. In addition these cuts are quite correlated with the cuts on 6-quark rapidities. By the 
dotted horizontal and vertical lines we have also marked limitations of the detectors planned 
at 220 and 420 meters from the collision point by the ATLAS and CMS collaborations. One 
can see that a pair of the 220 m detectors is not sufficient to measure Higgs boson. Both 
detectors on both sides are needed to measure most of the yield. 



D. Other backgrounds 

Our analysis in the previous subsection has been concentrated on the irreducible back- 
ground only. Other contributions, although, in principle, reducible, can in practice be also 
rather troublesome [l3[. The cross section for the gluonic dijets was found to be much larger 
than that for the bb jets j9[. If the gluon jets are misidentified as fe-jets, which was estimated 



to be 1.3% for the ATLAS detector, they contribute to the Higgs background. In Ref. [13 
the authors discuss in addition pile up events when the measured protons are not related 
to the exclusive Higgs production. Table 2 in their analysis presents detailed results for 



the issue. Inclusive double-pomeron processes [53| can also contribute to the background. 
Further analyses, especially for the Standard Model Higgs boson production, seem to be 
necessary to understand whether the Higgs boson can be identified in the exclusive produc- 
tion, perhaps not only in the bb decay channel. The present parton level analysis should 
be supplemented in the future by an additional analysis of bb jets by including a model of 
hadronization. Then standard jet algorithms could be imposed and the quality of the b and 
b kinematical reconstruction could be studied in detail. 



E. Some other remarks 

We have not been interested here in the precise estimation of the cross section but rather 
in understanding the signal-to-background ratio which is of the major importance for the 
upcoming Higgs boson searches in exclusive mode at the LHC. Consequently, we have pre- 
sented results with only one UGDF. This ratio is practically the same for other UGDFs. 
The absorption effects have been included here in a simple multiplicative form. They are 
expected to be the same both for the signal and the background, and thus are not affecting 
the ratio under consideration. The same gap survival factor has been used in both cases. 

As was mentioned above, in the current analysis we do not take into account the next- 
to-leading order QCD corrections in hard subprocess parts in both the bb background and 
Higgs CEP. Calculations of such corrections in the hard subprocess g*g* — > qq within the 
fc^-factorization approach are rather cumbersome, and we postpone them for our future 
studies. 

We have already analyzed the sensitivity of the results on the choice of UGDF. Different 
PDFs used to calculate UGDFs are defined in different range of factorization scales (gluon 
transverse momenta squared), some like CTEQ and MRST only for higher scales (q± min > 
1 GeV 2 ), some like GRV and GJR for lower values (q\ yTnin > 0.4 GeV 2 ). 
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Let us analyze how important are the low gluon transverse momenta in evaluation of the 
cross section. In Fig. [32] we show how the total and differential in M bb cross sections depend 
on the lowest value of the screening gluon transverse momentum squared used in evaluating 
the corresponding amplitude. There is much stronger dependence of the background than of 
the signal. This is caused by the specificity of matrix elements and the different three-body 
and four-body kinematics. It is interesting to note that at high lowest limit (> 1 GeV 2 ) the 
cross section for different gluon distributions coincide. This shows that the differences of 
the cross section between different UGDFs come mainly from the region of relatively small 
values of the screening gluon transverse momenta. There is a stronger sensitivity on q\ min 
for larger values of bb invariant mass (see the right panel in Fig. 132]) . 




FIG. 32: Total bb and Higgs CEP cross section as a function of the lowest cut on the screening 
gluon transverse momentum scale (left panel) and the invariant mass distribution of the bb CEP 
for different values of the cut and different factorization scales (right panel) . 

In Fig. [33] we show the invariant mass distributions of exclusive bb pair production for 
\ b Xl = ++ and H — (anti)quark helicity contributions with the realistic lower cut p b ± > 
40 GeV (so the high-p^ limit is concerned) on both quark and antiquark jets transverse 
momenta. The H — contribution clearly dominates, however the ++ contribution is not 
negligible, especially for very large invariant masses of the bb pair. Decreasing the cut-off on 
Pb± relatively enlarges the ++ contribution, making it important for low-p_i_ jets production. 



VII. CONCLUSIONS 

We have derived leading-order formula for the amplitude for EDD production of heavy 
quarks in the fcj_-factorization approach. This formula takes into account both gluon virtu- 
alities (transverse momenta) as well as the quark masses neglected in earlier works in the 
literature. We have shown that corresponding g*g* — y qq vertex is gauge invariant. We have 
also discussed purely QED double-photon component. 

Using the 2 — y 4 diffractive amplitude we have calculated differential cross section for 
cc and bb central exclusive production (CEP) in (anti)quark rapidities, quark and proton 
transverse momenta, transverse momentum of the qq pair and in azimuthal angles between 
outgoing protons and quark dijets in the whole four-body phase space for the nominal 
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FIG. 33: Invariant mass distributions of EDD bb pair production for different (anti)quark helicities 
with the lower cut (pi,j_ > 40 GeV) on both quark and antiquark jets transverse momenta. 



LHC energy y/s = 14 TeV. Large cross sections have been found in contrast to previous 
expectations in the literature. 

We have also discussed how the cross sections depends on quark masses. While at low 
quark-antiquark invariant masses the cross section for light quarks (u, d, s) is considerably 
larger than for heavy quarks (c, b) at large invariant masses the situation reverses. For 
instance, at invariant mass ~120 GeV (relevant for Higgs searches) it is the bb contribution 
which dominates. Since experimentally one can misidentify the other non-6 quark jets as 
6-jets, our calculation shows that this is not so dangerous provided that the misidentification 
probability is not too high. The gluonic jets seems in this context more difficult because of 
much larger cross section 13 . 

We have also calculated differential distributions for exclusive Higgs production as well 
as for b and b quarks (antiquarks) from the decay of the Higgs boson. We have used, for 
the first time in exclusive Higgs case, the vertex function which is consistently with the k±- 
factorization approach i.e. takes into account the gluon virtualities in the hard subprocess 
vertex. We have discussed the role of the off-shell effects. In contrast to the exclusive Xc 
production, the off-shell effects for Higgs boson are rather small and can be neglected given 
other sizeable theoretical uncertainties. 

The bb EDD and QED continua constitute an irreducible background to the exclu- 
sive Higgs boson production. We have discussed in detail how to improve the signal-to- 
background ratio by imposing cuts in quark rapidities, proton transverse momenta, longi- 
tudinal momentum fraction of outgoing protons. The analysis in the (yb, yi)-space is very 
useful to separate the two contributions as there they are located in quite different parts of 
this space. An optimal two-dimensional cut was proposed and the corresponding invariant 
mass distribution of the signal and background was presented. 
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